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Abstract. Echoes arise when external manipulations to a system induce a reversal 
of its time evolution that leads to a more or less perfect recovery of the initial state. 
We discuss the accuracy with which a cloud of trajectories returns to the initial state in 
classical dynamical systems that are exposed to additive noise and small differences in 
the equations of motion for forward and backward evolution. The cases of integrable 
and chaotic motion and small or large noise are studied in some detail and many 
different dynamical laws are identified. Experimental tests in 2-d flows that show 
chaotic advection are proposed. 
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1. Introduction 

Echoes arise when through suitable manipulations in a system the dynamics is reversed 
and a more or less complete recovery of the initial state is achieved. Acoustical echoes 
arise from reflections of sound at walls, spin echoes from reversals of magnetic fields 
(Hahn 1950, Carr and Purcell 1954), current echoes through a sequence of suitable 
electromagnetic pulses (Niggemeier et al 1993) and Loschmidt echoes from a reversal 
of momenta in a Hamiltonian system (Loschmidt 1876). That echoes can also appear 
in many particle systems is at first surprising since it seems to be in conflict with the 
irreversibility implied by the second law of thermodynamics. Closer inspection shows, 
however, that the recovery of the initial state is not perfect, and studies of the deviations 
tell a lot about the mechanisms that break reversibility. 

Several aspects of echo phenomena in dynamical systems have recently been studied 
in connection with Loschmidt echoes in quantum systems (triggered by Pastawski et 
al 1995 and Levstein et al 1998). An initial state |0) is propagated forward in time 
with Hamiltonian H and then back in time with a slightly different Hamiltonian H' (as 
suggested by Peres 1984). The loss of coherence is measured in terms of the fidelity 
^ e iH't/n e -iHt/h\^Y This is the same as calculating the overlap between the states 
\t) = e~' lHt / h \0) and \t') = e~ lHt '^ n \0), obtained by propagating the same initial state 
|0) for the same time t but with two different Hamiltonians. The decay of the overlap 
as a function of time and difference in Hamiltonian and the various time regimes have 
been the subject of several recent papers, e.g. (Jalabert and Pastawski 2001, Jacquod, 
Silvestrov and Beenakker 2001, Tomsovic and Cerruti 2002, Benenti and Casati 2002, 
Wisniacki and Cohen 2001, Prosen and Znidaric 2002, Prosen and Seligman 2002). 

The present paper is devoted to Loschmidt echoes in classical dynamical systems. 
If there is no difference between forward and backward equations of motion, the initial 
state is recovered perfectly. But what happens if there are small differences or if the 
system is exposed to noise? And what are the differences between echo experiments in 
integrable and chaotic systems? These questions will be addressed for Gaussian densities 
in linearized flows: they provide a convenient and sufficiently general class of densities 
in which a large variety of dynamical behaviour can be identified. 

Besides the obvious connection to the quantum echo experiments, the calculations 
are of some relevance for two other, directly classical situations: numerical trajectory 
reversals and reversibilty in advection. 

When the equation of motion are reversed in numerical calculations trajectories will 
typically not return to their starting point. For chaotic systems, the inherent sensitivity 
to initial conditions suggest an exponentially large deviation. A lack of growth has 
been used as an indicator for quantum regularity (Casati et al 1986). Obviously, no 
such problems should arise for perfect reversals and perfect numerical integrators since 
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the solutions to the equations of motion are uniquely specified by the initial conditions 
(baring singular points in the differential equations). The fact that trajectories do not 
return to their starting points thus reflects numerical inaccuracies from finite time steps 
and limited resolution. Using Gaussians to represent a cloud of initial conditions, noise 
to reflect truncation errors and differences between forward and backward integration 
algorithm this numerical reversibility experiment can be connected to the problem 
considered here. 

A popular and impressive demonstration of echoes in classical systems is provided 
by flow reversals in viscous liquids: a blob of dye can be stretched out until it is barely 
visible but upon reversal of the flow it reforms almost completely! The multimedia fluid 
mechanics CD (Homsey et al 2001) contains several demonstrations in laminar flows. 
The connection to chaos comes through experiments on chaotic advection (Aref 1984, 
2002, Ottino 1990). For instance, in their experiments on chaotic advection, Chaiken 
et al (1986) noted that if the dye passed through a chaotic region the recovery was less 
perfect, but they did not investigate this in detail. Since many of the typical flows can 
be realized experimentally, also the dynamics of the Gaussians discussed below should 
be experimentally accessible. 

The types of system considered here are classical. They may be exposed to additive 
white noise and forward and backward motion may differ. The state of the system is 
characterized by a smooth density in phase space, and the evolution equation is the 
Fokker-Planck equation with appropriate drift term. The discussion will be limited 
to densities that are Gaussian in shape, for which a reasonably complete general 
discussion is possible. Superpositions of such Gaussians can be used to approximate 
other densities. General expressions in arbitrary dimensions will be given, but their 
implications are most evident in 2-d conservative flows: this is also the case that is 
accessible in hydro dynamical systems. 

The outline of the paper is as follows. In the next section the dynamics of the 
center of mass and the variances for Gaussian densities will be discussed. In section 3 
this information is applied to the discussion of echoes in integrable systems (section 3a) 
and chaotic ones (3b). Section 4 contains a discussion of the results and some remarks 
on experimental tests in 2-d advection systems. 

2. Gaussian densities 

2.1. Outline of echo experiments 

The calculation of an echo naturally divides into two steps: the forward evolution up 
to some time T under one set of equations, followed by the backward evolution under 
a perhaps slightly different set of equations for the same time interval. If a mapping of 
initial conditions under a class of time evolutions can be found, say Pf — U\(T)pi for 
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the forward evolution under flow u 1; then we can write for the backward evolution with 
flow u 2 the formal expression p b = £/ 2 (— T)pf, so that the mapping to the echo state p e 
becomes 

p e = U 2 (-T)U 1 (T) Pi . (1) 

Thus, on the technical level it suffices to find U for the forward dynamics of a sufficiently 
large class of densities and systems. Note that in the presence of noise or in a dissipative 
system there is a difference between (i) comparing the initial state with the state 
obtained by propagating an initial condition over the complete cycle of forward and 
backward evolution, and (ii) comparing the states obtained by propagating the same 
initial condition forward in time with two different flows: dissipation and noise break 
the reversibility that permitted the change in protocol in the quantum case. 

The evolution equation for the densities p(x, t) is the Fokker-Planck equation, 

p = -V(up) + DAp (2) 

where D is the molecular diffusion constant. Echoes are most easily identified when 
initial and final density are sufficiently similar, as is the case for strongly localized 
objects. More complicated initial densities may be approximated by superpositions of 
localized ones. The dynamics for localized densities splits, in leading order in moments, 
into two parts, the motion of the center of mass and the changes in shape and size, as 
measured by the variances. This expansion may be extended to higher order moments 
of the density, but the equations become too cumbersome to analyze. 

2.2. Center of mass motion 

The center of mass of a localized density follows a classical trajectory xp(t), where 

xp(f) = u(x P (t),t). (3) 

The notation here is borrowed from hydrodynamic advection (Aref 1984, 2002), where 
u is a velocity field and x P the trajectory of a particle advected by the fluid. In other 
situations the velocity field u(x, t) has to be replaced by the right hand side f of an 
evolution equation x = f(x, t). 

For the backward propagation, where a small modification of the flow field is 
permitted, the trajectory may differ from the one during forward evolution and we 
need to estimate the differences between the two. Let xp(t) be a trajectory in a velocity 
field u and xp(t) + q(t) one in the perturbed velocity field u + <5u. To first order in q 
the equation becomes 



q(t)=A(t)q(t)+5u(xp(t),t), 



(4) 



where A is the linearization of the full velocity field u + Su at the trajectory x P , 

With the help of the monodromy matrix M(t), the solution to 

M = AM (6) 
with initial condition M(0) = 1, a formal solution can be given, 

q(i) = M(i) (g(0) + jf* rfrM- 1 (r)5u(x P (r), r)) . (7) 

This is the general solution for the displacement in a perturbed velocity field. A 
discussion of specific examples will be deferred to section 3 below. 

2. 3. Variances 

With x = x — xp(£) the coordinates relative to the center of mass, the density can be 
written as p(x, t) = p(x — xp(t),t) and the Fokker-Planck equation becomes 

p(x, t) = -V(u(x P + x, t)p) + (u(xp, t) • V)p + DAp , (8) 

where now all spatial derivatives are with respect to the relative coordinate x. The 
localization of the densities allows a linearization of the velocity field near the trajectory, 
i.e. 

u(x,i) = u(x P (i),f) + A(t)x (9) 

with the derivative matrix A, eq. (5). Then 

p = —(tiA)p — ((5c r A T ) ■ V)p + DAp . (10) 

For conservative Hamiltonian systems and incompressible flows, the trace of A vanishes; 
the discussion will henceforth be limited to that case. The next step in the analysis is 
to note that the equations are second order with a linear position dependence at most, 
so that a solution in terms of Gaussian densities is possible (Eckhardt 1990). In an 
n-dimensional phase space with variance matrix T(t) they are given by 

p(x, t) = 7T-"/ 2 (det IT 1/2 e-* Tr ~ lic (11) 

where T satisfies 

r = 2D + AT + TA T (12) 

(note that T as a kernel for a quadratic form is symmetric). With the help of the 
monodromy matrix M(t) also a closed expression for the variance matrix can be given, 

F(t) = M(t) (Vi + 2D J* (M(t) t M(t)) ^ dr^j M(t) T . (13) 

This is the central formula for the dynamics of the variances on which the calculation 
of the various situations can be based. For the specific cases studied in the next section 
a direct solution of (12) was found to be simpler. 
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3. Echoes 

3. 1 . Classical fidelity 

The analysis of echoes now proceeds as follows: We start with an initial density pi with 
variance matrix ry This density evolves under the influence of a velocity field u and 
additive white noise for some time T. At the end of this time interval the position of 
the density is xp(T) and the variance is Tf(T), as given by (13). Then the field is 
reversed. If the reversal is perfect, the new velocity field u' equals — u and the center of 
mass returns exactly. If there is a slight deviation the center of mass will move along 
a trajectory with a small displacement q(t) according to (7). In both cases, however, 
the variance does not return exactly, unless all noise is suppressed, i.e. D — 0. This 
holds quite generally, the only approximation being the Gaussian shape of the density 
and linerization of the velocity field near the trajectory. 

A general Gaussian echo may then be displaced with its center of mass by q and will 
have a variance matrix T e . For a comparison between initial and echo density of states 
the positions and variances can be used directly. But it is also possible to mimick the 
quantum fidelity expression and to introduce an overlap between classical phase space 
densities. The main difference is that quantum wave functions are normalized within 
the L 2 norm and classical densities are not. The proper definition of the fidelity as the 
cosine of the angle between the two densities in Hilbert space then is 



= /rfx g 2 (x)Pe(x) m . (14) 

(Jdxp i (x)2) 1 / 2 (J dxpe(x )2 ) V2 



Without the normalization the overlap could change even if both densities evolve in the 
same way. The definition given by Prosen and Znidaric (2002) thus has to be modified. 
For the case of two Gaussians, an initial density pi centered at zero and with variances 
Tj, and an echo density 

p e = TT-/ 2 (det r e )" 1/2 e -(x-«0 T r.- 1 (x-q) (15) 
the overlap becomes 



O = 2 n/2 t 



Vdet iydet T± T( r -i +r -i) q 

\ det(r, + r e ) • 



The overlap integral indicates two very different kinds of contributions: The 
prefactor measures the reduction in overlap by spreading of the density. The exponential 
factor accounts for the rapid drop off in overlap when the centers of the Gaussians are 
separated; the Gaussian form is clearly connected to the Gaussian tails in the density, 
and would be different, e.g., for exponential tails in the density. 
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3.2. Shear flows 

In laminar flows neighboring trajectories see slightly different velocities and separate 
linearly in time. When combined with noise a cubic growth of the variance results. 
Specifically, consider a 2-d shear flow 

«=(?) (17) 
of shear rate a. The associated monodromy matrix is 

"<«>=(;?)■ («) 

During forward evolution the variances become 

rj{ } = rg + 2aT§T + a 2 r^T 2 + 2L>T + -a 2 DT 3 (19) 
Ti? = rg + aT§T + aDT 2 (20) 

rg ) = rg + 2DT. (21) 

The T 3 contribution to the variances has also been discussed by Rhines and Young 
(1983) in the context of fluid mixing. 

For the reversal we allow for a different shear rate a' and some perturbation in the 
velocity field. If q/ is an initial displacement in the trajectory and 5u = (ui,u 2 ) a 
constant perturbation to the velocity field, the displacement of the echo will be 

Qi,e = qi,f - («i - <x'q 2 ,f)t - a'u 2 t 2 /2 (22) 
<?2, e = Q2j ~ OLU 2 t - a'u 2 t 2 /2 . (23) 

The variances become 

rff = rg + 2(a - a')r$T + (a- a'fT^T 2 

+ 4DT + J (2a 2 + 8a' 2 - 6aa')DT 3 (24) 

rg = rg + (a - aOrgT + (a - 3a')^T 2 (25) 

rg = rg + 4/JT. (26) 

This result for the variances may be verified for a few limiting cases: (i) Without shear 
a = a' = the diagonal elements increase like 4DT as for regular diffusion over a 
time interval 2T. (ii) Without diffusion (D=0) the determinant of the matrix does not 
change, (iii) Without diffusion (D=0) and equal shear in the forward and backward 
direction (a = a') the reversal is perfect and the initial variances are restored, (iv) The 
parameters a' = —a correspond to the situation that the backward integration is just 
a continuation of the forward integration and the expressions (24)-(26) agree with (19)- 
(21) for an evolution time of 2T. (iv) For equal shear a = a' in forward and backward 
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direction the variances are 

= rg + ADT + -a 2 DT 3 (27) 
3 

rg = rg - 2aDT 2 (28) 

rg = rg + 4L>T; (29) 

the prefactor of the cubic term in rff is smaller than would be obtained from (24) for 
a time 2T, indicating that the reversal of the shear induced broadening is partial and 
not complete. 

The different time regimes in the variances are easily identified. Consider the terms 
linear in time first: diffusion will be noticable on a time scale Tjj ~ 1/-D, the differences 
in the shear rates on a time scale T$ ~ l/\a — a'\. Thus, a large diffusion can swamp 
the effects from the difference between the two Hamiltonians. For the nonlinear terms, 
the one with the difference in shear rates appears around T§, and the one with the cubic 
term in diffusion near 1/a. In typical applications the latter should be smaller than T$. 

The classical fidelity contains the determinants of the variances. In the absence of 
diffusion, D = 0, we have det T e = det Tj and 

det(r; + r e ) = 4 det T t + (a - a') 2 T 2 T$ (30) 

so that the fidelity has a 1 — const • T 2 behaviour for short times and a l/(\a — a'\T) 
decay for times longer than T$. With diffusion and equal forward and backward shear, 
a = a', we have 

det T e = det Ti + 4(rg + T^)DT 

+ IQD 2 T 2 + 4T^aDT 2 + -T ( $a 2 DT 3 + -a 2 D 2 T 4 (31) 

3 3 

and 

det (Ti + r e ) = 4detr l + 8(rg + rg)/xr 

8 4 
+ 16D 2 T 2 + SrflaDT 2 + -T^^DT 3 + -a 2 D 2 T 4 (32) 

so that a short-time behaviour 1 — const • t and a long-time behaviour of l/(\faDT) 
follow. 

The contributions from displacements enter in the exponent and can in principle 
introduce rapid decays. Consider, e.g., the case a = a' and weak diffusion, weak shear 
and short times, so that DT and aT are smaller than the initial variances. Then 
p( e ) ~ p(») anc l tn e only time dependence will come from the perturbations qo and 5u: 
the exponent will contain a polynominal of fourth order in time, and this leads to a 
rapid and faster than exponential decay. If diffusion is added the increase in variance 
can compensate part of the displacement growth, but for strong diffusion and large times 
an exponential decay will remain: the square of q 2 will increase like T 4 and the variance 
increases like T 3 , so that the ratio increases linearly, giving an exponential decay. 
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Thus, without displacement the overlap integral decays algebraically as contained 
in the prefactor, but with displacement the decrease can be dramatic when the densities 
are separated by more than their widths. 

3.3. Chaotic systems 

As a model for a chaotic system take a simple hyperbolic motion, u = (Ax, —\y) and 
assume an initial density aligned with the unstable (x-) and stable (y-) direction, i.e. 
= 0. After forward evolution we have 

r£? = r«e 2AT + ^(e 2A '-l) (33) 
A 

rg^rge-^ + ^l-e- 2 ^). (34) 

Thus, there is an exponential contraction down to the limit set by diffusion. If the 
backward integration has a slightly different stretching rate A', then for the echo 

r (e) = r g e 2(A-A')T + ^ (e 2AT _ l)& -2WT + _ g _ 2VT) ^ 

rg = rge- 2 ( A - A ') T + j (1 - e- 2AT )e 2A ' T + ^(e 2A ' T - 1) . (36) 

As in the previous case we can study various limiting situations, such as A, A' — > 0, 
where linear diffusion results, or A = —A', where the expressions (33) and (34) for times 
up to 2T are recovered. When forward and backward stretching rates are the same, 
A = A', the variances become 

rS e 1 ) = rg+2^(l-e" 2AT ) (37) 
A 

r 2 2 ) = rg+2^(e 2A ' r -l). (38) 

Note that the variance in x has hardly changed whereas the one in y grows exponentially. 
The reason is that the x-variance grows during the forward integration and collapses 
then during the backward evolution, down to a limit set by the diffusional broadening. 
The exponential growth and contraction is thus almost perfectly compensated. For 
the ^/-variance we have first the contraction, down to the limit set by diffusion. The 
expansion during the backward evolution then starts from this finite amplitude, and 
not from the exponentially small contraction of the deterministic evolution of the initial 
variance. As a result, the echo is stretched out along the direction that was the stable 
one during forward evolution. If forward and backward evolution are interchanged, then 
so is the orientation of the spreading of the density: it will then point in the x-direction. 
For this growth to be noticable the time evolution has to be followed for times longer 
than about (ln(Ar«/D))/(2A). 
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The behaviour of the classical fidelity for short times is quadratic or linear, for 
D = and D ^ 0, respectively. On longer times there is an exponential decay, like 
exp(— |A — A'|T) without diffusion and like exp(— X'T) with diffusion. 

The appearance of differences between Lyapunov exponents reflects a relation 
between forward and backward flow: the stable and unstable manifolds are aligned. 
A more general difference between forward and backward flow will break this alignment 
and introduce exponentials in A'. 

Small perturbations in position and in the velocity fields will grow exponentially. 
Eq. (7) gives 



The contributions from the displacement to the decay are weaker than in the linear shear 
flow, since the stretching of the variances is in the same directions as the separation of 
trajectories, so that an exponential increase in q 2 can be compensated by an exponential 
increase in variances. However, in the absence of diffusion and for the same Lyapunov 
exponents in forward and backward direction an exponentially growing displacement can 
lead to a drastic drop off, like exp(— expAt), in overlap, simply because the Gaussians 
are shifted relative to each other. 

4. Final remarks 

Already the simple examples in the previous section show a wide range of dynamical 
behaviour in classical echo experiments. There are two ingredients: the variation in 
variances and a displacement between initial and echo density. Without displacement 
the decay in the integrable system is slower than in the chaotic one, and the same applies 
when there is a displacement in both systems. However, the overlap can drop off faster 
in an integrable system with displacement compared to a chaotic one without. 

The drop off from the displacement is connected with the Gaussian shape of the 
densities: if the tails fall off more slowly then also the overlap will decay more slowly 
as a function of displacement. The shape dependence should be stronger in integrable 
systems than in chaotic ones, since the stretching of variances and shapes goes in parallel 
with the exponential growth of displacement. 

The calculations are based on Gaussian densities and linearizations of the flow fields 
near trajectories. This becomes questionable if the densities spread out too far, a 
problem that occurs more likely and more quickly in chaotic systems than in integrable 
systems. In integrable systems the dangerous terms come from quadratic dependencies 
of the winding frequencies on action and from curvatures introduced when mapping the 
tori onto position space. The density will then coil up in whirls. In a 2-d of freedom 




(39) 
(40) 



<?2,e = (q 2 ,f + U 2 / \')e 



u 2 /X' ■ 
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system this can be expected to happen linearly in time. In chaotic systems the linear 
approximation is applicable if the stable and unstable manifolds are close to straight lines 
within the area covered by the density and if the variations in stretching and contraction 
rate are small. Typically, this will limit the time interval to a logarithmically short one. 
For longer times the density develops tendrils that follow the manifolds as they wiggle 
through phase space (for an early discussion of such effects within the quantum maps, 
see Berry et al, 1979). In both cases the degree to which the Gaussian approximation 
breaks down depends on details of the nonlinear contributions and has to be considered 
for specific models. 

The transformation to a comoving frame eliminates the center of mass motion and 
emphasizes the linearized dynamics near the trajectory. The same discussion thus 
applies to stationary points. Two examples, shear flow or parabolic type and hyperbolic 
type, have been discussed here, the third class, elliptic type, has periodically oscillating 
variances. 

Two-dimensional flows provide an ideal testing ground for the results presented 
here. Localized spots of dye can be prepared as initial densities and their motion can 
be followed in various 2-d flows (Aref 1984, 2002, Ottino 1990, Homsey 2001, Chaiken 
1986). Elliptic and hyperbolic points can be realized most easily in cellular flows (Jiitner 
et al 1997, Williams et al 1997, Rothstein et al 1999) Molecular diffusivities are fixed 
by the selection of dye and solvent, but variations of shear rate through the amplitude 
of the velocity fields and initial variances through the size of the spot provide enough 
degrees of freedom to explore the full range in behaviour. In particular the 2-d Lorentz 
force driven flows (Jiitner et al 1997, Williams et al 1997, Rothstein et al 1999) should 
have enough flexibility to study echoes in flows and to provide quantitative tests of 
the various expressions derived here. It should be possible to see the cubic laws in the 
variances in integrable systems, the chaos assisted spreading in chaotic systems and the 
dependence on the order of forward/backward propagation. 

It would also be of interest to reconsider the experimental protocol of Chaiken et al 
(1986): they followed a closed line of dye. After reversal the line was displaced left and 
right of the original trace and also had slightly different widths. This suggests different 
histories of the center of mass of small line elements and passage through regions with 
different degrees of chaos. Using the complete tracings of the flow field as in (Voth, 
Haller and Gollub 2002) it should be possible to pin down the different regions and 
interactions and to characterize the reversal completely. 
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